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A  REVISED  FORMULA  FOR  THE  CALCULATION  OF  GAMMA -RAY 


INTRODUCTION 


In  earlier  publications  an  empirical  formula  was  developed  for  predicting 
the  effectiveness  of  two-legged  air  ducts  through  concrete  for  the  attenuation 
of  gamma  radiation.1 < ^  Results  were  directly  applicable  to.  the  case  of 
radiation  streaming  through  entranceways  and  air  ducts  leading  into  protective 
structures. 

During  the  past  two  years,  however,  certain  events  have  taken  place 
which  make  it  desirable  to  revise  the  earlier  formula.  Chapman  has  modified 
and  improved  his  computer  technique  for  calculating  the  streaming  of  gamma 
radiation  through  ducts.3  Also,  Chapman  has  made  additional  measurements  on 
the  streaming  of  gamma  rays  from  Au1’®,  Cs13?,  and  Co® 

In  view  of  the  above  developments,  a  reappraisal  was  made  of  the  earlier 
empirical  formula  in  the  light  of  information  presently  available.  This  note 
is  concerned  with  the  development  of  a  new  empirical  formula.  Some  sample 
problems  are  worked  to  demonstrate  its  usefulness.  The  accuracy  and  limita¬ 
tions  of  the  formula  are  discussed.  A  statistical  analysis  is  used  to  show 
that  the  formula  has  essentially  no  bias  and  a  quite  small  random  error. 


DEVELOPMENT  OF  THE  FORMULA 


The  formula  as  originally  developed  to  approximate  attenuation  factors 
in  two-legged  concrete  ducts  was1** 


D  $ 

£  -  0.25 


0.907 


.2.864 


.  27534.  2.667  0.710 
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where 

D0  =  dose  rate  in  mr/hr  at  a  distance  from  the  source  of  one  foot  in  air 

D  =  dose  rate  in  mr/hr  inside  the  shelter 

H  »  height  of  entranceway  in  feet 

W  =  width  of  entranceway  in  feet 

Lj>  length  in  feet  of  first  leg  of  duct 

L2*  length  in  feet  of  second  leg  of  duct 

EQ*  average  energy  of  source  gamma  radiation  in  Mev 

The  formula  was  considered  to  be  valid  when  the  following  inequalities 
are  true 

0.662  SE0«  3.000  Mev 
1.0  £  H  *  6.0  feet 
1.0  £  W  £  6.0  feet 
2  £  Lj  s  36  feet 
1  £  H/W  £  2 
Lj/H  *  6 
L2/H  *  6 
lq/W  i  2 
L2/W  *  2 

The  formula  had  a  safety  factor  already  incorporated  into  it,  ao  that 
one  could  be  95-percent  certain  that  the  actual  attenuation  factor,  D/0o,  is 
at  least  aa  small  as  the  attenuation  factor  calculated  by  the  formula. 

Because  of  the  availability  of  both  new  experimental  data  and  an  improved 
computerised  calculational  technique,  it  was  decided  to  seek  an  improved 
empirical  formula.  The  first  attempt  was  to  fit  data  to  a  formula  of  the  type 


K  (B  ♦  AH)  *(W  ♦  flp 


(2) 
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where  K  and  the  a4 a  are  adjustable  parameters.  The  width  and  height  were 
incremented  by  an  amount  depending  on  EQ,  to  account  for  the  fact  that 
scattering  of  gamma  rays  does  not  occur  at  the  surface  but  at  some  depth 
below  the  surface  of  the  scatterer.  This  technique  of  mathematically 
incrementing  the  physical  dimensions  of  the  duct  was  successfully  used 
by  Chapman  in  his  computer  calculations.^ 

It  was  soon  found,  however,  that  Equation  2  did  not  give  an  appreciably 
better  fit  than 


"l  “2 
(H/W)  W 


(3) 


Since  Equation  3  is  simpler  than  Equation  2,  it  was  decided  to  forego 
the  refinement  of  incremented  cross  sectional  dipensions.  Curve  fitting 
yielded  the  following  formula 


2  ■  0.244 
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.  2.701  2.705  _  0.581 
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(4) 


Comparison  of  the  above  with  Equation  1  shows  minor  changes  in  the  adjust- 
sble  parameters.  These  changes  result  from  addition* 1  information  which 
has  become  available  since  Equation  1  waa  postulated. 

Study  of  Equatim  4  indicated  further  simplification  might  be  in  order. 
The  W  in  H/W was  multiplied  out  and  other  rounding  off  of  exponents  was 
performed,  producing 
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It  la  noted  that  Equation  5  ia  not  dlmenalonally  correct,  alnce  the  left 
hand  aide  la  dimenalonlsaa  while  the  right  hand  aide  has  units  feet”2** 
Mav“0*®.  This  failing  ia  a  common  fault  of  empirical  equations.  Although 
tnnaerically  correct  results  are  obtained,  there  remalna  the  conceptual 
difficulty  of  unequal  units.  One  should,  then,  imagine  that  the  right  hand 
side  of  Equation  5  la  multiplied  by  a  constant  K  which  has  magnitude  unity 
and  dtmenslone  of  feet2**  Msv° • c .  The  fact  that  such  units  are  not  phys¬ 
ically  meaningful  ia  a  consequence  of  the  fact  that  the  equation  is  an 
empirical  one  which  has  not  been  derived  from  purely  physical  cons ids rat Iona, 


Since  a  statistical  analysis  indicated  that  Equation  5  fit  the  available 
data  almost  as  well  as  Equation  4  fit  the  data,  Equation  5  is  adopted,  and  it 
is  recommended  for  calculations  of  attenuation  of  gamma-ray  dose  in  two-legged 
air  ducts  through  concrete  whenever  the  following  inequalities  hold: 

0.3  3.7  Mev 

1.0SVSHS6.0  feet 

2.0  SL-  i  24.0  feet,  i  ■  1,  2 

Li/H  S6.  ,  i  ■  1,  2 

1.0  £  H/W  £  2.0 

It  will  be  noted  that  the  inequalities  listed  above  are  somewhat  differ¬ 
ent  from  the  conditions  listed  earlier  for  the  validity  of  Equation  1.  It 
is  believed  that  the  set  of  inequalities  specified  for  Equation  5  are  more 
reasonable  for  the  case  of  protection  from  fallout  radiation. 

In  later  sections,  the  curve-fitting  procedure  will  be  described  and 
the  accuracy  of  Equation  5  will  be  discussed. 


CURVE  PITTING 

Having  decked  to  fit  an  equation  of  the  form  of  Equation  3,  the  next 
step  was  to  select  the  data  to  be  fit  and  to  perform  the  curve-fitting 
analysis. 

Fifty  experimental  measurements,  corresponding  approximately  to  those 
reported  by  Chapman, ^  were  chosen.  Also,  Chapman's  computer  code  was  used 
to  calculate  the  same  cases.  This  gave  an  additional  56  datum  points, 
because  two  computer  calculations  were  performed  for  each  of  the  six  cases 
involving  the  gamma  rays  of  Ha21*,  since  Na24  has  two  well  separated  gamma-ray 
energies.  In  addition,  117  fictitious  cases  were  selected  at  random  from  the 
domain  of  interest.  Those  cases  were  calculated  using  Chapman's  code  to 
provide  another  117  datum  points.  Since  agreament  between  Chapman's  calcula¬ 
tions  and  experiment  is  excellent  in  most  cases,  it  seemed  reasonable  to 
believe  that  his  calculations  could  be  trusted  to  be  accurate  anywhere  within 
the  domain  of  energy  and  duct  geometry  where  it  has  been  tested. 

The  method  tor  selection  of  the  117  fictitious  cases  consisted  of  choos¬ 
ing  lengths,  widths,  and  heights  at  random  from  the  domain  of  interest;  i.e., 
subject  to  the  stipulated  inequalities.  Thus  all  allowable  dimensions  were 
equally  likely  to  be  chosen. 

Since  there  is  little  or  no  experimental  data  for  energies  above  3.0 
Mev,  and  since  energies  resulting  from  fallout  are  in  the  neighborhood  of 
1  to  2  Msv,  it  wes  decided  to  limit  the  energy  range  to  from  0.3  Mav  to  3.7 
Mev,  with  concentration  around  2.0  Mev.  This  was  accomplished  by  choosing 
energies  at  random  from  a  population  of  energies  with  a  maximum  density  at 
2.0  Mev  and  with  a  density*  falling  oft  linearly  from  the  maximum  to  aero  at 
0.3  Mev  and  3.7  Msv.  A  histogram  of  the  energies  actually  selected  is  shown 
in  Figure  1. 


Using  the  223  cases  described  as  data,  the  IBM-1620  computer  program, 
"Stepwise  Regression  Analysis  Program"  (STRAP)  Revised*  was  used  to  perform 
the  analysis.  The  parameters  were  found  to  be 


K  « 

“1  * 
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The  formula  then  became 


0.244  a3  «  2.701 
0.9758  an  •  2.705 
3.048  a5  =  0.581 
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(4) 


The  223  cases  were  then  recalculated  using  Equation  4  and  the  results  compared 
with  the  measured  results  and  those  calculated  with  Chapman's  code. 

Figure  2  shows  a  histogram  of  the  ratios  obtained  when  the  223  results  of 
Equation  4  were  divided  by  the  223  datum  points.  The  mean  was  found  to  be 
1.011,  and  the  standard  deviation  was  found  to  be  0.128.  Also  on  Figure  2  is 
a  theoretical  Normal  density  function  having  those  values  for  mean  and 
standard  deviation. 

A  x2  goodness -of -fit  test  was  performed.  The  value  of  x2  Obtained  was 
6.91.  A  x3  6.4  would  disprove  the  assumption  of  normality  at  the  50-per¬ 
cent  confidence  level,  while  »  x  »  9.0  would  disprove  normality  at  the 
75-percent  confidence  level.  Thus,  the  x2  analysis  provided  no  reason  to 
disbelieve  the  assumption  of  normality.  Therefore,  normal  theory  was  used  to 
make  predictions  regarding  the  population  of  all  possible  ducts  within  the 
domain  of  interest. 

Next,  the  ratioa  between  datum  points  and  formula  results  were  obtained 
for  the  simpler  formula 
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Figure  3  shows  the  histogram  of  the  ratioa  obtained.  While  the  mean, 
1.007,  more  oloaely  approximated  the  assumed  mean  of  1.000,  the  standard 
deviation  wee  somewhat  larger,  0.133.  Thua,  approximately  95  percent  of 
the  cases  lie  withln*27  percent  of  the  Chapman  solution  and/or  the  measured 
results,  and  approximately  two-thirds  are  within  ±13.5%. 

Since  Equation  5  ia  simple  and  gives  satisfactory  agreement  with  the 
data,  Equation  5  ia  recommended  for  calculations  of  gamma-ray  streaming 
through  concrete  ducts. 


c; 
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USING  THE  FORMULA 

In  order  to  obtain  conservative  (safe)  estimates  of  the  dose  rate 
inside  a  shelter,  use  is  made  of  the  formula 
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2.7  „  0.6 
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A  simple  way  is  now  described  for  using  the  formula. 

Enter  Figure  4  with  the  value  of  Lj  for  the  duct  of  interest,  and  find 
the  corresponding  value  for  LJ  *  .  Similarly,  find  L2^*^,  also  from  Figure  4. 

Find  E00*6  from  Figure  5.  If  the  energy,  E0,  is  unknown,  assume 
E0  -  2  Mev. 

Now  substitute  into  Equation  6  to  find  a  conservative  estimate  of  the 
attenuation  factor  of  the  entranceway. 

It  should  be  noted  that  all  the  discussion  up  to  this  point  has  been 
concerned  with  a  point  isotropic  gamma-ray  source  located  on  the  centerline 
axis  just  at  the  outside  opening  of  the  entranceway.  In  a  radiation  shelter, 
however,  what  is  really  of  interest  is  the  protection  factor;  i.e.,  the  ratio 
of  the  dose  rate  at  3  feet  above  the  ground  in  an  unprotected  place  to  the 
dose  rate  in  the  shelter.  It  is,  therefore,  necessary  to  modify  Equation  6, 
in  order  to  obtain  a  conservative  estimate  of  the  desired  quantity,  PF, 
protection  factor.  Define 


Fg  • 


(7) 


where  D/D0  is  given  by  Equation  (6). 

The  protection  factor  is  then  defined  as 


PF  ■ 


1 

PgFS 


(8) 


where  Fa  is  a  source  factor.  For  fallout  radiation,  reference  is  made  to 
Figure  6  (due  to  J.  C.  LeDoux).  When  the  first  leg  is  horiaontal 


For  a  vertical  first  leg,  where  the  aperture  covering  ia  contaminated, 

,  2  . 

P*  *  h  Ah 
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For  a  vertical  first  leg,  where  the  aperture  covering  ia  not 
contaminated, 


Fs  Nl.  A 

1  a 


The  three  different  cases  are  shown  in  Figure  6,  which  also  giv«3 

values  for  the  factors  A  ,  A.  ,  and  A  . 

v  n  a 

SAMPLE  PR0B1EMS 

The  simplified  empirical  formula  will  now  be  used  to  ralculate  dose 
rates  within  three  sample  ducts. 

Case  1 

Consider  a  1-foot  square  duct  with  Lj  *  3.5  feet  and  L2  *  2.0.  The 
source  energy  is  1.25  Mev,  corresponding  to  Co^O. 

Equation  5  ia  now  used: 


)  ,  _ (1)  (1) 

>o  4(3.5  *  2.0)2'7(1.25)C 


■  0.114  x  10“ 


The  above  can  be  obtained  by  computer  calculation,  by  slide  rule,  or  by 
logarithms. 

For  the  sane  case,  Chapman’s  computer  code  gives  the  answer: 


g  •  0.1223  x  10"2 
o 


Tamil  measured  this  duct  situation  as: 


E  •  0.108  x  10"2 

uo 


Case  2 

*  i*  ’ 

Consider  a  6-foot  -square  duet  with  L}  •  12  feet  and  L2  ■  IS  feet. 
The  source  energy  la  again  1.25  Mev, 

Equation  5  gives 


D  _  6*6 

K  *  . .  ~-e-w—  ■  -  w  *  m  0.385  x  10 

Do  4(12  •  15)2*r(1.25)07® 
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Chapman’s  code  gives 


D  -4 

g  «  0.4068  *  10 

o 

The  measured  value  was? 


£  >  0.354  x  10~** 


Case  3 

As  a  final  case,  consider  a  3-foot  square  duct  with  Lj  ■  L2  ■  7.5  feet. 
The  source  energy  is  0.662  Mev,  corresponding  to  Cs^?, 

Equation  5  gives 


2  ■  0.163  x  10"3 


Chapman's  code  gives 

£  *  0.1662  x  10"3 
o 

The  measured  value  waa4 

£  ■  0.156  x  10"3 


Thus  It  ia  seen  that  In  «U  three  cates  there  it  good  agreement  between 
the  simple  formula  and  both  Chapman' s  computer  calculations  and  experimental 
values. 


CONCLUSIONS 


It  has  been  shown  that  the  attenuation  factor  of  a  two-legged  air  duct 
through  concrete  can  be  satisfactorily  fit  by  the  empirical  equation 
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for  fallout  gamma  radiation  over  the  domain  of  duct  dimensions  of  interest. 

The  formula  is  a  good  one  in  the  sense  that  there  is  no  apparent  bias 
in  its  predictions.  It  does,  however,  have  c  standard  deviation  of  13.5 
percent.  Therefore,  a  safety  factor  should  be  incorporated  to  assume  that 
predictions  of  dose  rates  inside  shelters  will  be  unlikely  to  be  low.  A 
recommended  conservative  formula  is 
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Using  the  above  formula,  one  can  be  95-percent  certain  that  the  actual 
dose  rste  will  not  be  greater  than  the  formula  predicts. 

In  order  to  obtain  the  protection  factor  in  a  fallout  field  use  is  made 
of  the  equation 


PP 
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The  above  equation  will  giva  a  conaervative  estimate  of  the  protection 
factor  of  a  fallout  shelter,  assuming  radiation  can  enter  only  through  the 
antranceway . 
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Figure  1.  Histogram  of  117  rendoofly  selected  energies. 


Figure  3.  Histogram  of  ratios  between 
formula  values,  Rjf  and  data  points,  Rq,  t or  Equation  5 
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